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Pion photo- and electroproduction has been studied at threshold and in the resonance region 
below W < 2 GeV. At threshold tt° production can be very well explained within a dynamical 
' model derived from an effective chiral Lagrangian. The final state interaction is nearly saturated by 

single charge exchange rescattering. In the resonance region new electroproduction data at Q 2 = 1 
GeV 2 has been analyzed with MAID and longitudinal and transverse photon helicity amplitudes 
' have been determined for different resonances. A detailed study of the E/M and S/M ratios of the 

■ N — > A transition shows a zero crossing of Rem near Q 2 = 4 GeV 2 , whereas the Rsm becomes 

CN , increasingly negative at large Q 2 . 

1 Introduction 

The unitary isobar model MAID is a model for single pion photo- and electroproduction off 
protons and neutrons [1]. It is based on a non-resonant background described by Born terms and 
vector meson exchange contributions and nucleon resonance excitations modeled by Breit-Wigner 
functions 

r-j . 71" — 1^ ' I - * ■ \ / 

Both parts, background and resonance are separately unitarized. This has been achieved by a 
K-matrix unitarization in the case of the background 



O 
=5 



gfiMAID) = exp (i5 a ) cos 6 a v% a (W, Q 2 ) , (2) 



^ , and by introducing a unitary phase cj)R for the resonance excitations 

t 7W [W,Q ) - A a {Q ) M 2 r _ W 2_ iMrTr e ■ W 

The phases 8 a are the elastic pion-nucleon scattering phases in a particular channel a = {l,j,t} 
below the inelastic threshold of two-pion production. In order to take account of inelastic effects, 
the factor exp (iS a ) cos 5 a is replaced by \{n a exp (2i5 a ) + 1] with the inelasticity parameters i] a at 
higher energies. Additional background terms are included to account for S- and P-wave pion loop 
effects. 

In the case of the Dynamical Model (DMT) [2], the background contribution is given by 



& a (DMT) = e i&a cos<5 Q 
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with the full ttN scattering reaction matrix R^- In this case the pion loop effects that are especially 
important near threshold are generated dynamically and show up as a principal value integral over 
the reaction matrix. 
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Presently, we have included 8 nucleon resonances, 3 Deltas: ^33(1232), S < 3i(1620) and 1)33(1700) 
and 5 N*s: Pn(1440), L>i 3 (1520), S"n(1535), Sn(1650) and Fi 5 (1680). All of them are included 
with longitudinal and transverse electromagnetic couplings. The corresponding helicity amplitudes 
A 1 / 2 , A 3 / 2 and Si/ 2 can be fitted to experimental data and can be freely changed in the MAID 
program. 

We will show results obtained with the Dynamical Model for the threshold region and present 
some recent fits to newer electroproduction data with MAID for the resonance region with a 
discussion of the Q 2 evolution of the E/M and S/M ratios of the N — > A(1232) transition. 

2 Results in the threshold region 

For 7T° photoproduction, we first calculate the multipole Eq + near threshold by solving the coupled 
channels equation within a basis with physical pion and nucleon masses. The coupled channels 
equation leads to the following expression for the pion photoproduction t-matrix in the TT°p channel: 

V>(W) = v ln o(W) +v ln0 (W)g 7T o p (W)t w o p ^ 7T o p (W) 

+ v ln+ (W)g 7T+n (W)t w+n ^ 7T o p (W), (5) 

where t^Op^^Op and t 7T + n ^ 7T o p are the irN t-matrices in the elastic and charge exchange channels, 
respectively. They are obtained by solving the coupled channels equation for irN scattering using 
the meson-exchange model of Ref. [3]. Our results for ReEo + show that practically all of the final 
state interaction effects originate from the n + n channel and mainly stems from the principal value 
integral of Eq. (5). In this approach the t n N matrix contains the effect of irN rescattering to all 
orders. However, we have indeed found that only the first order rescattering contribution, i.e. the 
one-loop diagram, is important. This indicates that the one-loop calculation in ChPT is a reliable 
approximation for 7r° production in the threshold region. 

If the FSI effects are evaluated with the assumption of isospin symmetry (IS), i.e., with averaged 
masses in the free pion-nucleon propagator, the energy dependence in Re Eq + in the threshold region 
is very smooth. Below ir + threshold the strong energy dependence (cusp effect) [4] only appears 
because of the pion mass difference and, as we have seen above, is related to the coupling with 
the 7r + n channel. In most calculations, the effects from the pion mass difference below the tt + 
production channel are taken into account by using the K-matrix approach [5], 



Re EH = Re E ol ( IS ) ~ a *N lo c Re E%$ (IS) J 1 - — , (6) 

where to and oo c are the ir° and ir + cm. energies corresponding to W = E p + w 7 and m n + m n +, 

respectively, and a n N = 0.124/771,,-+ is the pion charge exchange threshold amplitude. ' (IS) 
is the 7r°' + photoproduction amplitude obtained with the assumption of isospin symmetry (IS), 
i.e., without the pion mass difference in Eq. (4). Such an approximation is often used in the data 
analysis in order to parametrize the Eq + multipole below 7r + n threshold in the form of Eq + (E) = 
a + by/l — (uj/uj c ) 2 . Numerically this approximation is very precise and differs only by 10% at the 
7T° threshold and becomes indistinguishable above the ir + threshold. In Fig. 1 the results obtained 
within this approximation scheme are represented by the solid curve and compared to the ChPT 
calculation (dash-dotted curve) [6]. Over the whole energy range the difference is rather small and 
within the experimental uncertainties. Huge effects, however, arise if the cusp would be neglected 
or obviously if the FSI effects would be totally ignored (dotted curve). 

In Fig. 1 we also compare the predictions of our model for the differential cross section with re- 
cent photoproduction data from Mainz [8,9]. The dotted and solid curves are obtained without and 
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Figure 1: Real and imaginary parts of the Eq + multipole and differential cross sections below and 
above ir + threshold for 7p — ► 7r°p. The solid and dashed curves obtained with and without the 
cusp effect, respectively. The dotted curve is without FSI. The dash-dotted curve is the result of 
ChPT [5]. Data points for £0+ from Mainz [8] (A), [9](o) and Saskatoon [10](»), and for da from 
Mainz [8](«), [9](o). 



with FSI effects, respectively. It is seen that both off-shell pion rescattering and cusp effect substan- 
tially improve the agreement with the data. This indicates that our model gives reliable predictions 
also for the threshold behaviour of the P- waves without any additional arbitrary parameters. 

Pion electroproduction provides us with information on the Q 2 = —k 2 dependence of the 
transverse Eq+ and longitudinal Lq + multipoles in the threshold region. The "cusp" effects in the 
Lq+ multipole is taken into account in a similar way as in the case of Eq + , 



ReLH = ReLH (IS) - a wN co c Re (IS) «/l - ^ , (7) 

V W c 

where all the multipoles are functions of total cm. energy W and virtual photon four-momentum 
squared Q 2 . It is known that at threshold, the Q 2 dependence is given mainly by the Born plus vector 
meson contributions in v^ n: as described in Ref. [1]. Similar to pion photoproduction, the K-matrix 
approximation and full calculation agree with each other within a few percent. In Fig. 2 we show 
our results for the cusp and FSI effects in the Eq + and Lq + multipoles for 7r° electroproduction 
at Q 2 = 0.1 (GeV/c) 2 , along with the results of the multipole analysis from NIKHEF [11] and 
Mainz [12]. Note that results of both groups were obtained using the P-wave predictions given by 
ChPT. However, there exist substantial differences between the P— wave predictions of ChPT and 
our model at finite Q 2 . To understand the consequence of these differences, we have made a new 
analysis of the Mainz data [12] for the differential cross sections, using our DMT prediction for 
the P-wave multipoles instead. The 5-wave multipoles extracted this way are also shown in Fig. 2 
by solid circles. We see that the results of such a new analysis gives Eq + multipoles closer to the 
NIKHEF data and in better agreement with our dynamical model prediction. However, the results 
of our new analysis for the longitudinal Lo+ multipoles stay practically unchanged from the values 
found in the previous analyses. Note that the dynamical model prediction for Lq + again agrees 
much better with the NIKHEF data. Further details are given in Ref. [13] 

In contrast to DMT, in MAID the FSI effects are taken into account using the K-matrix 
approximation, namely without the inclusion of off-shell pion rescattering contributions (principal 
value integral) in Eq. (4). As a result, the S-, P-, D- and F- waves of the background contributions 
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Figure 2: Real parts of Eq + and Lq + for ep — ► e'ir°p at Q 2 =0.1 (GeV/c) 2 . Notations are the same 
as in Fig. 1. Data points from NIKHEF [11](°) and Mainz [12](A). The results of the present work 
obtained by using the P- waves of our model are given by (•). 



are defined as 

ta(MAID) = exp (iS a ) cos 6 a v^(q E , k). (8) 

However, as we have found above, dynamical model calculations show that pion off-shell rescatter- 
ing is very important at low pion energies. The prediction of MAID for Eq + {tt°p) at threshold, 
represented by the dotted curves in Fig. 1 lies substantially below the data. It turns out that it 
is possible to improve MAID, in the case of tt° production at low energies, by introducing a phe- 
nomenological term and including the cusp effect of Eq. (6). In this extended version of MAI D2000, 
we write the Eq + (tt°p) multipole as 

Re E^l° = Re (MAID98) + E cusp (W, Q 2 ) + E corr (W, Q 2 ) , (9) 

where 



E cusp (W,Q 2 ) = -a wN iu c ReE^l {MAID98) J 1 - — . (10) 



uj- 



The phenomenological term E corr which emulates the pion off-shell rescattering corrections (or 
pion-loop contribution in ChPT) can be parameterized in the form 

E corr (W,Q 2 ) = (1 + ^ V)2 ^(Q 2 )' ( n ) 

where Fjj is the standard nucleon dipole form factor. The parameters A and B are obtained by 
fitting to the low energy ir° photoproduction data: A = 2.01 x 10 _3 /m 7r + and B = 0.71 fm. 

3 Results in the resonance region 

For pion photoproduction in the resonance region we have recently performed a fit of MAID for 
both the low-energy A region and the medium-energy resonance region up to W = 1700 MeV, 
where resonance parameters have been obtained as a part of the BRAG 1 partial wave benchmark 
analysis [14]. In Fig. 3 we present a new fit of preliminary results on electroproduction, p(e, e'p)ir® , 

1 Baryon Resonance Analysis Group, http://cnr2.kent.edu/~manley/BRAG.html 
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measured by the JLab Hall A collaboration [15]. The data has been taken at backward angles at 
Q 2 = 1.0 GeV 2 in the cm. energy range from 0.95 GeV to 2.0 GeV. With a dataset of 363 data 
points in 3 observables, da = dax + edai^daiT and daxr and pion angles of 146,151 and 167 
degrees we performed a data analysis with MAID. 

From the 20 possible resonance parameters we have varied 18 by fixing E 2 ~ and S2- of the 
.Di3(1520) because we did not find enough sensitivity in the data which are only taken at backward 
angles. In table 1 we give the result of our fit both for the multipole and the helicity amplitudes. 
The multipole amplitudes are compared to the default values of MAID. For the A(1232) resonance 
we give in addition the E/M and the S/M ratios. Both are consistent with the previous MAID 
fits to photo- and electroproduction [16]. The Rsm ratio is very well determined by the daLT data 
and shows the tendency to larger negative values for increasing Q 2 , while the Rem ratio is much 
more uncertain and also the model uncertainties are larger than for the S/M ratio. From daiT 
we also find a large sensitivity to the So+ amplitude of the 511(1535) resonance in the minimum 
around W = 1500 MeV as well as for the S2- amplitude of the 1)33(1700) resonance in the second 
maximum around W = 1650 MeV. Furthermore most of the structure in da and in cIxttt above 
W = 1700 MeV is explained by the Mi- amplitude of the Z?33(1700) resonance. However, the 
MAID model does not include higher resonances so far. 
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Figure 3: Preliminary experimental results from the JLab Hall A collaboration [15] at Q 2 = 1.0 
GeV 2 , 6 n = 167° and e = 0.9 as functions of the cm. energy W. The dashed lines show the standard 
MAID2000 calculations and the solid curves are the results of the fit to the data. 

In order to get information on resonance properties a careful analysis has to be taken. First of 
all, a partial wave decomposition in terms of multipoles is necessary to get information on quantum 
numbers as angular momentum, spin, parity and isospin. Second, a background separation is needed, 
as especially in pion production a large background is produced by the strong pion nucleon coupling. 
However, nucleon Born terms and vector meson exchange contributions are not the only source of 
background. Loop effects can give very large contributions especially for S- and P-waves, the most 
famous example is the Eq + in tt° photoproduction at threshold. Our resonance extraction is based 
on the imaginary parts of the full multipoles in a specific spin-isospin channel at the resonance 
position. This minimizes the model dependence since resonance positions are mostly well known 
and this method can be applied to any given partial wave analysis. It is very similar to the method 
that has been applied for the A(1232), however, as Watson's theorem is no longer fulfilled for 
higher resonances, some uncertainties have to be accepted. The helicity amplitudes extracted by 
this way are "dressed" amplitudes and contain contributions from vertex corrections. They could 
be undressed by subtracting background contributions and unitarization corrections. By such a 
procedure they could be directly related to the Breit-Wigner couplings A^iQ 2 ) of Eq. (3), however, 
such a method will always be model dependent. 

The helicity amplitudes Ai/ 2 ,^3/2 and Si/ 2 are determined from the pion electroproduction 
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multipoles at the resonance position 
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The equivalent photon energy k^- and the pion momentum are given in the cm. frame and 
evaluated at the resonance position, where also the pion electroproduction multipoles are obtained, 
A = ImA(W = Mr) for A = E, M, S. For the transverse amplitudes these formulas agree with 
PDG and Ref. [17]. For longitudinal amplitudes we found different definitions in the literature, here 
we use a definition consistent with notations used in DIS [18]. 
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Figure 4: The Q 2 dependence of the magnetic G* M form factor and the E/M and S/M ratios 
at W = 1232 MeV. The solid and dashed curves are the MAID and dynamical model results, 
respectively. JLab results discussed here are shown at Q 2 = 1 (GeV/c) 2 , Bates results at 0.126 
(squares), the Mainz double polarization result of S/M at 0.12 (full circle), data at 2.8 and 4.0 
from Ref. [19] (squares) and preliminary results of the ratios in the range of 0.1 - 0.8 are from 
Bonn [20] (open circles). Results of our analysis for the ratios at 2.8 and 4.0 are obtained using 
MAID (full circles) and the dynamical models (triangles). In the case of G* M they fully agree with 
Ref. [19]. For older data of G* M see Ref. [16]. The photoproduction results of Mainz [21] are placed 
at the lowest Q 2 value. All numbers are given in units of (GeV/c) 2 . 

In Fig. 4 we show our extracted values for the magnetic form factor G* m /3Gd and the ratios 
Rem and Rsm together with other data determined in different ways in recent experiments and 
data analyses on a semidog scale. For photoproduction the Mainz results [21] are shown, around 
0.125 GeV 2 the values of the Bates analysis and the result of the Mainz measurement with recoil 
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MAID 
defaults 


MAID 
fit results 


helicity 
amplitudes 


P 33 (1232) E^ 2 

53/2 
°1+ 

Rem 
Rsm 


-0.44 
20.1 
-1.31 
-2.2 
-6.5 


-0.38± 0.17 
19.0± 0.4 

-1.55± 0.05 
-2.0± 0.9 
-8.1± 0.2 


A 1/2 -72± 3 
A s/2 -135± 3 
51/2 18± 1 


Pi i (1440) M\' 2 

~1 /o 

5^1 


1.87 
1.26 


2.2± 0.2 
0.6± 0.1 


Ai /2 -59± 9 
5i /2 -11± 4 


Di 3 (1520) 

~ 1 /o 

5^ 


-0.05 
1.71 



-0.05 
1.2± 0.2 



A 1/2 -47± 6 
A a/2 26± 4 
5i /2 


=r T7o 

511(1535) So 1 ; 2 

c.l/2 


4.37 
0.1 


3.5± 0.8 
-1.2± 0.1 


A 1/2 61± 14 
5i /2 -15± 2 


5 31 (1620) El' 2 

C°/ 
°0+ 


0.21 
0.89 


3.4± 0.3 
1.7± 0.5 


A 1/2 -50±5 
5i /2 -18± 7 


5n(1650) El' 2 

c.l/2 


2.8 



3.0± 0.4 
-1.2± 0.6 


A 1/2 43± 6 
5i /2 -12± 6 


Pis (1680) E L 3 '_ 2 
Ml' 2 

S£ 


0.32 
0.83 



-0.06± 0.03 
0.80± 0.05 
-0.10± 0.03 


A 1/2 -52± 9 
A 3/2 33± 9 
5i /2 -7± 2 


D 33 (1700) P 2 i/2 
Ml' 2 
S 3 ' 2 


-0.85 
0.30 



-1.0± 0.2 
1.1± 0.1 
0.2± 0.1 


A 1/2 104±12 
A 3/2 -4±12 
5i /2 -14± 7 


PV-PS mixing: A m 


450 


350 ± 35 





Table 1: Proton resonance multipoles (A = lmA(W = M r ) in 10~ 3 / m 7r), helicity amplitudes (in 
1(T 3 GeV~ 1/2 ) and values of the PV-PS mixing parameter A m (in MeV) as in MAID2000 and 
obtained in our fit at Q 2 = 1.0 GeV 2 . The _Di3(1520) E2- and §2- amplitudes were fixed. The 
E/M and S/M ratios of the A (1232) are given in percentage. 

polarization, in the medium Q 2 range the preliminary data of Bonn [20] and at high Q 2 our analysis 
of the JLab Hall C data [19] with MAID and the DMT. The main difference between our results 
and those of Ref. [19] is that our values of Rem show a clear tendency to cross zero and change 
sign as Q 2 increases. This is in contrast with the results obtained in the original analysis [19] of the 
data which concluded that Rem would stay negative and tend toward more negative values with 
increasing Q 2 . Furthermore, we find that the absolute value of Rsm is strongly increasing. 

4 Summary 

With the unitary isobar model MAID and the dynamical model DMT we have two very good 
tools available to analyze data of pion photo- and electroproduction and to plan new experiments 
with increased sensitivity to specific questions. While DMT includes pion loop contributions that 
are especially important for low partial waves (S and P), MAID originally was constructed only 
from tree diagrams, resonance excitations and K-matrix unitarization contributions. Therefore, to 
get better agreement for S-waves, MAID has been extended phenomenologically by low-energy 
corrections and the unitary cusp effect. At higher energies the PS-PV mixing of MAID that was 
already introduced from the beginning, also serves for this purpose to effectively taking into account 
of loop contributions. 

With all parameters already fixed by ttN scattering and pion photoproduction in the resonance 
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region, DMT describes pion photo-production at threshold very well, similar to the calculations in 
ChPT or with dispersion relations. For electroproduction at threshold we also find good agreement 
with the experiment at Q 2 = 0.1 GeV 2 , however, we also have problems describing the recent 
Mainz data at Q 2 = 0.05 GeV 2 as it also appears in ChPT calculations. In the resonance region 
both models DMT and MAID can equally well describe photo- and electroproduction data up 
to W = 1700 MeV by fitting the photon helicity amplitudes A1/2, A 3 / 2 an d £1/2 of the individual 
resonances. Here, we have demonstrated this with the preliminary Hall A data of JLab at Q 2 = 1.0 
GeV 2 . Even with a dataset limited to backward pion angles we are able to determine quite a 
few resonance parameters in satisfactory precision and can also give longitudinal couplings where 
previous information practically did not exist. For the E/M and S/M ratios of the Delta resonance 
we combine our previous fits of Mainz, Bates and JLab Hall C data with our new analysis and find 
a consistent Q 2 evolution of these ratios with a slowly rising Rem that crosses zero around Q 2 = 4 
GeV 2 , and for the longitudinal coupling a Rsm that significantly increases to larger negative values 
at high Q 2 . If expectations from pQCD will be fulfilled, a sharp rise in the E/M ratio towards 
100% should be seen in the next generation of experiments above Q 2 = 5 GeV 2 and a leveling of 
the S/M ratio to a constant value. 
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